home *** CD-ROM | disk | FTP | other *** search
/ io Programmo 60 / IOPROG_60.ISO / soft / c++ / gsl-1.1.1-setup.exe / {app} / src / msvc / demo / main.c < prev   
Encoding:
C/C++ Source or Header  |  2002-04-18  |  1.0 KB  |  47 lines
  1. #include <stdio.h>
  2. #include <gsl/gsl_rng.h>
  3. #include <gsl/gsl_vector.h>
  4. #include <gsl/gsl_blas.h>
  5.  
  6. /* See the file README.txt for information on compiling this program */
  7.  
  8. int
  9. main (void)
  10. {
  11.     const gsl_rng_type * T = gsl_rng_default;
  12.  
  13.     gsl_rng * r = gsl_rng_alloc (T);
  14.  
  15.     int i;
  16.  
  17.     printf ("Here are ten random numbers in the range 0-99:\n");
  18.  
  19.     for (i = 0; i < 10; i++)
  20.     {
  21.         int k = gsl_rng_uniform_int (r, 100);
  22.         printf(" %d", k);
  23.     }
  24.  
  25.     printf("\n");
  26.  
  27.     {
  28.         double x[5] = { 1.0, 2.0, 3.0, 4.0, 5.0} ;
  29.         double y[5] = { 5.5, 4.4, 3.3, 2.2, 1.1} ; 
  30.         gsl_vector_view v = gsl_vector_view_array(x, 5);
  31.         gsl_vector_view w = gsl_vector_view_array(y, 5);
  32.  
  33.         printf("blas operation DAXPY\n");
  34.  
  35.         printf("x:"); for (i = 0; i < 5; i++) { printf(" %g", x[i]); } ; 
  36.         printf("  y:"); for (i = 0; i < 5; i++) { printf(" %g", y[i]); } ;
  37.         
  38.         gsl_blas_daxpy (1.0, &v.vector, &w.vector);
  39.         
  40.         printf("  a x + y:"); for (i = 0; i < 5; i++) { printf(" %g", y[i]); } ;
  41.     }
  42.  
  43.     printf("\n");
  44.  
  45.     return 0;
  46. }
  47.